AD-A151  414  EFFECTS  OF  FINITE  CURRENT  CHANNEL  WIDTH  ON  THE  CURRENT  1/ 
CONVECTIVE  INSTRBILITVCU)  NAVAL  RESEARCH  LAB  WASHINGTON 
DC  P  SATVANARAVANA  ET  AL.  BS  MAR  85  NRL-HR-5494 
UNCLASSIFIED  F/G  28/9  NL 


AD-A151  414 


-  * 

••• 


Effects  of  Finite  Current  Channel  Width 
on  the  Current  Convective  Instability 

i 

P.  Satyanarayana  and  J.  Chen 

•  Science  Applications,  Inc. 

McLean,  VA  22102 

M.  3.  Keskinen 

Geophysical  and  Plasma  Dynamics  Branch 

Plasma  Physics  Division  v 


March  5, 198S 


This  research  was  partially  sponsored  by  the  01 
under  Subtask  S99QMXBC,  work  unit  001{ 


of  Naval  Research  and  the  Defense  Nuclear  Agency 
nd  work  unit  title  *Plasma  Structure  Evolution.* 


4 


NAVAL  RESEARCH  LABORATORY 
Washington,  D.C. 


Approved  for  paMie  re  lew;  distribution  unttmitod. 

85  08  07 


DT1C 

ELECTE 


.aSS-  CA'ON  --  S  '^Gc 


REPORT  DOCUMENTATION  PAGE 


i  3£?Q*~  SECwRlTV  C.ASS  c'CAT  ON 

UNCLASSIFIED 


O  RESTRICTIVE  MARK'NGS 


:<*  sscuRitv  class. f.cation  autK)R»tv 


:d  DECLASS-F'CATlON.  DOWNGRADING  schedule 


j  Distribution,  ava.lability  op  report 

Approved  for  public  release;  distribution  unlimited. 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

NRL  Memorandum  Report  5494 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a  NAME  CP  PERFORMING  ORGANIZATION 

Naval  Research  Laboratory 


6o  office  symbol 

( if  applicable) 

Code  4780 


7a  NAME  OF  MONITORING  ORGANIZATION 


**  ADDRESS  <Oty,  Stare,  and  Zip  Code) 

Washington,  DC  20375-5000 


7b  aDORESS  (Gfy,  Stare,  and  ZIP  Code) 


8a.  NAME  OF  PUNOiNG  /  SPONSORING 
ORGANIZATION 

ONR  and  DNA 


8b  OFFICE  SYMBOL 
(If  applicable) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


Re  ADDRESS  (City,  State,  and  ZIP  Code) 

Arlington,  VA  22217  Washington,  DC  20305 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT  J 

ELEMENT  NO 

NO 

NO 

ACCESSION  NO 

(See  page  ii) 

1  title  (include  Security  Classification} 

Effects  of  Finite  Current  Channel  Width  on  the  Current  Convective  Instability 


V  PERSONAL  AUTHOR(S) 

Satyanarayana,  P.  *  Chen,  J  *  and  Keskinen,  M.J. 


3a  TYP£  OP  REPORT 

Interim 


1 3b  TIME  COVEREO 

from  10/83  to  10/84 


14  DATE  OF  REPORT  (Tear,  Month,  Day) 

1985  March  5 _ 


S  PAGE  COUNT 

41 


6  SUP°LEMENTAR''  NOTATION 

♦Science  Applications,  Inc.,  McLean,  VA  22102 


r?  COSATI  CODES  I 

f-eld 

GROUP 

SUB-GROUP  ‘  • 

(Continues) 


18  SUBJECT  rERMS  [Continue  on  reverse  if  necessary  and  identify  by  block  number ) 

Plasma  instability  Fluid  instability 

High  latitude  ionosphere  Finite  current  channel  width  ,  ^  - 


*9  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

The  effects  of  the  finite  current  channel  width  on  the  current 
convective  instability  are  studied  both  analytically  and  numerically. 
First  using  a  sharp-boundary  field-aligned  current  distribution  which  has  a 
finite  width  along  the  plasma  density  gradient,  the  dispersion  relation  is 
obtained  analytically.  It  Is  found  that,  for  the  long  wavelength  modes 
(kyd  «  1)  where  the  nonlocal  effects  are  most  prominent,  the  growth 

rate  y  Is  proportional  to  l^d/L2 1  in  the  collisional  limit  f  »  |cu  1 1  , 


(Continues) 


.’0  DiSTPlBUNON/ AVAILABILITY  OF  ABSTRACT 
CQuNCLASSIFlEO/UNLIMITEO  □  SAME  AS  RPT 


□  OTIC  USERS 


2  T  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


22.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

J.  D.  Huba 


22b  TELEPHONE  (Include  Area  Code) 

(202)  767-3630 


22c,  OFFICE  SYMBOL 

Code  4780 _ 


OO  FORM  1473,34  MAR 


3i  APR  edition  may  be  used  until  exhausted 
All  other  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


1 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO. 

NO. 

NO. 

ACCESSION  NO. 

62715H 

RR033-02-44 

DN880-024 

61153N 

DN230-291 

16.  SUPPLEMENTARY  NOTATION  (Continued) 

This  research  was  partially  sponsored  by  the  Office  of  Naval  Research  and  the  Defense  Nuclear  Agency  under 
Subtask  S99QMXBC,  work  unit  00102  and  work  unit  title  “Plasma  Structure  Evolution.” 

19.  ABSTRACT  (Continued) 

where  d  is  the  half-width  of  the  current  channel,  L  is  the  plasma  density 
gradient  scale  length,  Vd  is  the  field-aligned  current  velocity,  vin  is  the 
ion-neutral  collision  frequency  and  <u  is  the  perturbation  frequency.  For 
the  long  wavelength  modes  in  the  inertial  limit  fv  <<  |u|),  the  growth 
rate  y  scales  as  jv^d/ L  |  a,  where  a  *  1/2  (2/3)  for  k^/k ^  ^e/ve£  much  less 
than  (greater  than)  |<u|/n  and  k  (k  )  is  the  wavenumber  parallel 

i  z  y 

(perpendicular)  to  the  magnetic  field,  0  fll.)  is  the  electron  (ion) 

e  l  <  <  ■' 

gyrofrequency ,  and  v  is  the  electron-ion  collision  frequency Numerical 
results  are  also  presented  for  a  diffuse-boundary  current  velocity 
distribution.  Applications  to  the  high  latitude  ionosphere  are  discussed. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


ii 


CONTENTS 

I.  INTRODUCTION  . 

II.  MODE  STRUCTURE  EQUATION . 

III.  DISPERSION  RELATION . 

Collisional  Limit . 

Inertial  Limit  . 

IV.  NUMERICAL  RESULTS  . 

V.  DISCUSSION  AND  SUMMARY . 

ACKNOWLEDGMENTS . 


REFERENCES 


EFFECTS  OF  FINITE  CURRENT  CHANNEL  WIDTH 
ON  THE  CURRENT  CONVECTIVE  INSTABILITY 

I.  INTRODUCTION 

From  a  variety  of  experimental  observations  [see  recent  reviews  by 
Fejer  and  Kelley,  1980,  Vickrey  and  Kelley,  1983,  and  Hanulse  et  al.  , 
1981],  it  is  now  known  that  the  high  latitude  ionosphere,  from  the  auroral 
zone  into  the  polar  cap,  is  a  highly  structured  and  nonequilibrium  medium 
containing  plasma  fluctuations  and  structures  with  scale  sizes  ranging  from 
hundreds  of  kilometers  to  meters.  Several  processes,  e.g.,  particle 
precipitation,  plasma  instabilities,  and  neutral  fluid  d3/namics,  have  been 
proposed  to  account  for  high  latitude  ionospheric  plasma  structure 
[Keskinen  and  Ossakow,  1983a].  The  plasma  instabilities  that  have  received 
Che  most  attention  thus  far  are  the  E_xB_  instability  [  Simon,  1963;  Linson 
and  Workman,  19  70;  Keskinen  and  Ossakow,  1982,  1983b;  Huba  et  al ;  1983, 

Keskinen ,  1984]  and  the  current  convective  instability  [ Lehnert ,  1958; 

Kadomtsev  arj  Nedospasov,  1960;  Ossakow  and  Chaturvedi,  1979;  Huba  and 
Ossakow,  1980;  Keskinen  et  al. ,  1980;  Chaturvedi  and  Ossakow,  1981,  1983; 
Satyanarayana  and  Ossakow,  1984;  Gary,  1984;  Huba ,  1984],  The  current 

convective  instability  results  from  the  coupling  of  a  magnetic  field- 
aligned  current  and  a  plasma  density  gradient  transverse  to  the  magnetic 
field  in  the  presence  of  electron-meutral  collisions  —  a  configuration 
likely  to  occur  in  the  auroral  zone  ionosphere  [Vickrey  et  al.  ,  1980].  All 
of  the  previous  studies  of  the  current  convective  instability  have  assumed 
that  the  width  of  the  current  transverse  to  the  magnetic  field  is 
infinite.  However,  recent  observational  results  [see,  for  example,  Bythrow 
et  al. ,  (1984)]  indicate  that  the  transverse  dimension  of  the  field-aligned 
current  distributions  may  be  comparable  to  or  less  than  the  density 
gradient  scale  length.  Since  the  instability  requires  both  a  density 
gradient  and  a  field-aligned  current,  the  properties  of  the  instability  may 
be  modified  if  one  uses  a  more  realistic  geometry  with  a  current  channel  of 
finite  width.  Indeed,  physical  intuition  would  lead  one  to  expect 
reduction  in  the  growth  rate  and  modification  in  the  stability  boundary  for 
narrow  current  channels.  In  this  analysis  we  remove  the  infinite-width 
current  simplification  and  study  the  effects  of  finite  current  channel 
width  on  the  current  convective  instabilty.  The  analytical  and  numerical 
Manuscript  approved  October  4,  1984. 
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results  will  show  that,  in  the  long  wavelength  limit,  k^d  <<  1,  where  the 
nonlocal  results  are  most  prominent,  the  growth  race  deoends  on  the  global 
quantities  such  as  ( V , d/L“),  where  V,  is  the  current  velocitv,  d  is  the 

half-width  of  the  current  channel  and  L  is  the  densitv  gradient  scale 

length.  These  intrinsically  nonlocal  effects  may  be  of  interest  for 

interpretation  of  the  observational  data.  We  discuss  the  nonlocal  results 
and  apply  them  to  the  high  latitude  ionosphere. 

The  organization  of  the  paper  is  as  follows.  In  Section  II  we  derive 
the  general  dispersion  relation  describing  the  mode  structure  of  the 
current  convective  instability  with  a  finite  current  channel  width.  In 

Section  III  we  present  a  theoretical  analysis  of  the  mode  structure 

equation  while  in  Section  IV  we  give  numerical  results.  Finally,  in 

Section  V  we  summarize  our  findings  and  discuss  applications  to  the  high 

latitude  ionosphere. 

II.  MODE  STRUCTURE  EQUATION 

In  this  paper  we  use  a  slab  model  shown  in  Fig.  (1).  The  plasma  is 

A 

assumed  inhomogeneous  along  the  north-south  direction  (x)  and  is  assumed  to 
support  equilibrium  currents  along  the  magnetic  field.  The  magnetic  field 
is  uniform  and  vertical  ( z )  with  the  field-aligned  currents  taken  to  be  of 

A 

the  form  J  =  J  (x)  z.  The  primary  objective  of  this  paper  is  to  confine 
—  o  — 

the  field-aligned  currents  to  a  finite  region  along  the  direction  of  the 

A 

density  gradient  (x)  and  to  study  the  effects  of  finite  current  width  on 
the  current  convective  instability.  We  neglect  weak  altitude  dependent 
density  gradients  as  well  as  east-west  (y)  gradients  and  transverse 
equilibrium  electric  fields.  The  temperature  effects  responsible  for 
diffusive  damping  are  ignored  as  is  magnetic  field  shear.  Electron  inertia 
is  neglected  since  current  convective  instability  is  a  low  frequency 
instability  ( m  «  ft  |.  Furthermore,  we  confine  ourselves  to  the  F-region  of 
the  ionosphere  and  ignore  the  electron-neutral  collision 
frequency  (ven)  compared  to  the  electron-ion  collision  frequency  fv^]. 
Due  to  the  variation  of  the  field-aligned  currents  and  the  background 
density  along  x,  we  perform  a  nonlocal  analysis  and  derive  a  mode  structure 
equation  for  the  perturbed  potential.  This  procedure  allows  us  to  study 
short  wavelength  and  long  wavelength  modes.  The  basic  two-fluid  equations 


describing  the  electron  and  ion  dynamics  in  the  rest  frame  of  the  neutrals 


are 


+  V 


(n  V  ) 

ct—  a 


0 


(1) 


V  x  B 

m  a  f—  +  V  •  V)  V  -  q  n  f  E  +  — - )  -  m  n  v  V  -  R 

a  a  3t  —a  — a  at  a  —  c  a  a  an-a  -a 


(2) 


where  a  denotes  species  (e  for  electrons;  i  for  ions),  n^  is  the  nlasma 

density,  _V  is  the  fluid  velocity,  is  the  mass  of  the  species,  qQ  is  the 

charge  of  the  species,  c  is  the  speed  of  light,  _E  is  the  electric 

field,  v  is  the  plasma  particle  collision  frequency  with  neutrals, 
an 

and  R  is  the  friction  force  given  by 


(3) 

(4) 

(5) 


V°.  -  O'*  (4) 

iz 

In  obtaining  these  equations,  terms  proportional  to  ven/ne  <<  1  have  been 

neglected.  By  setting  Soi  *  0*  it  can  be  easily  shown  that  «  0. 

We  choose  the  friction  force  to  be  Rg  *  -  so  that  momentum  conservation 

demands  v.  in  »  v  .m  and  thus  V?  -  0.  The  net  drift  velocity 
0  le  6  ei  e  iz  q 

V,  =  V  -  V.  then  becomes  V,  -  V  .  More  generally,  however,  neutral 

d  ez  iz  d  ez  o 

particles  also  participate  in  momentum  conservation,  and  V  need  not  be 
zero.  The  above  simplification  implies  that  the  electrons  are  the  primary 
current  carriers. 

In  order  to  derive  the  nonlocal  mode  structure  equation,  we  use  a 

A 

perturbation  of  the  form  f  -  f(x)  exp  [-ifmt-k  y-k  z)j,  where  u>  -  <d  +  iy , 

y  z  r 

Y  >  0  for  growth  and  k  and  k  are  the  mode  numbers  in  y  and  z 

y  z 

directions.  We  then  linearize  Eqs.  (1)  and  (2),  and  finally  we  subtract 
the  electron  continuity  equation  from  the  ion  continuity  equation  [Eq.  (1)] 
and  use  quasi-neutrality  condition  to  derive  the  mode  structure  equation. 
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We  consider  only  electrostatic  perturbations  and  do  not  perturb  the 
magnetic  field  [Ossakow  and  Chaturvedl,  1979].  By  perturbing  Eq.  (2)  for 
both  electrons  and  ions,  we  obtain  the  perturbed  velocities  in  the  z 
direction 


ez  m  v  . 

e  ei 


The  perpendicular  perturbed  velocities  are  given  as 

V  ,  -  x  2  -  Jv./f?  )E.1  (9) 

el  B  1  —  '  ei  e;  1J 

f  oi  +  iv  )  « 

lu-f  Hi  ‘ - 57 - III  (10) 

We  substitute  Eqs.  (7)-(10)  in  Eq.  (1)  and  use  _E  »  -  ^4  in  the 
electrostatic  limit.  Then  we  subtract  the  electron  continuity  equation 
from  the  ion  continuity  equation  and  impose  quasineutrality  to  obtain 

2*  * 

^-4  +  p(x)  +  q(x)<f>  -  0,  (ID 

3x  3X 

where 


p(x)  -  n;/nQ 


i^(Vvei]ni 

q(x)  *  -  k  - 

1“  +  lv-,n) 


kz  Vx) 


(<*>  *  kzvd(x))  +  ivm 


- -  Ik  fn'/n  }  -  i(n  /v  )  k^j 

■  1  y  0  0  '  p  p-f 


e  ei  zJ 


This  equation  agrees  with  the  nonlocal  mode  structure  equation  derived 
earlier  by  Satyanarayana  and  Ossakow  (1984)  and  Huha  (1984). 

The  effects  of  finite  current  width  are  contained  in  the  third  term  in 
Eq.  (13).  We  solve  Eq.  (11)  exactly  for  two  specific  profiles  of  the 
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current  velocity  V^(x).  In  the  next  section,  we  obtain  the  analytic 
dispersion  relation  for  a  "waterbag”  distribution. 


Vd,  Ixl  <  d 


0  ,  !x|  >  d 


where  2d  is  the  current  channel  width.  We  choose  an  exponential  plasma 

density  orofile  n  (x)  =*  n  exp  (x/L)  so  that  n"/n  *  1/L.  In  section  IV  we 
o  o  ^  oo 

present  the  numerical  results  for  both  the  waterbag  current  model  and  the 
diffuse  Gaussian  current  model. 


III.  DISPERSION  RELATION 


The  dispersion  relation  is  obtained  by  matching  the  logarithmic 
derivatives  of  the  wavefunction  inside  and  outside  the  waterbag  at  the 
boundaries,  x  =  +d,  and  requiring  that  the  solution  be  evanescent 

for  x  *  +  ®.  First,  we  will  cast  Eq.  (II)  into  a  Schrodinger  type 

equation  by  using  the  transformation 


■»  1 

<t  -  ^(x)  exp  [- j  J  p(x)dx] . 


r . A 


Then,  we  have 


where 


-Mf  +  Q2(x)  =  0 

3xz 


Q2(x)  =*  q(x)  -  d2/4  -  p"/2 


and  p'  *  3p/3x.  For  an  exponential  density  profile,  nQ  =  nQ  exp  (x/L), 
p  =  1/L  and  p'  =  0.  Thus,  for  the  waterbag  profile  V^(x)  given  by  Eq.  (14) 
and  for  the  exponential  density  profile,  Eq.  (16)  can  be  written  as 


d  vout  ~2  . 

^3—  -  °o  *out  -  0 

32^fn  2 

I?"  +  =  0 


|x|  >  d 


ix|  <  d 


A  * 

*  *  ^’  *.*  *.' 


with 


k~  +  l/f4L“l’ 
v  J 


i  k~'  Q  / v  .  'ft 
z  e  si  ■ 

fw  -*■  ivin' 


k2  Vd  gl 

fcu-k2Vd)  (ai  +  iuin) 


IVL  - 1  fW1' 


„  i  k"  f  Q  /v  ,  . 

kZ  +  1/(4L  1  + - 

■  y  (a.  +  iv.  ) 

V  i  r"* 


where  primes  indicate  derivatives  with  respect  to  x.  The  current 
convective  instability  is  a  nearly  purely  growing  mode  with 

;j  =  to  +  iv  -  iy.  This  will  be  shown  to  be  true  a  posteriori  in  the  next 

r 

section.  For  a>  <  y,  the  sign  convention  applicable  in  Eqs.  (18)  and 
(19)  is  such  that  Q2  >  0  and  Q2  >  0.  In  Eq.  (20)  the  sign  of  Vd/L  (  <  0)  is 

chosen  so  that  >  0  and  the  growth  rate  is  positive  for  positive  ky. 

Furthermore,  in  the  following  we  assume  that  I kz”7d  I  <<  I'111  in  the 
denominator  and  k2  (S^/v^)  <  k  /L.  First,  we  rewrite  ot2  and  Qc2  as 

follows: 


where 


k  f  ai  +  iv 


t“  +  ivin) 


k  (  Hi  -  0)+  ){ ID  -  ID_ 

“1“  +  iv<J 


=.  u+u_  =  fk2ky/k2)  |Vd/Lin1 
«  -  ( iv/2 ) l  1  +  (1  +  4 \<*lQh  ))  \ 


k  +  1/4 

y 
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and  k  ,  k  ,  and  d  are  the  normalized  quantities  k  L,  k,L,  and  d/L 
y  2  y  ** 

respectively.  We  note  that  ujj.  are  the  two  roots  of  the  local  dispersion 
relation  i  =  0  1.  Since  the  potential,  07(x) ,  in  Eq.  (16)  is  symmetric 
about  x  =  0,  we  can  restrict  our  analysis  to  eigenfunctions  of  definite 


parity. 


In  the  region  |x|  <  d  we  have 


A  cos  ?(Lx/L' 


for  even  parity 


=  B  sin  fO^x/L)  for  odd  parity 


where  A  and  B  are  constants.  In  the  region  |x|  >  d,  we  have  an 
exponentially  decreasing  solution  for  either  case 


ip  .  =  C  e 
Tout 


“  Qo  'I' 


where  C  is  a  constant.  The  matching  condition  in  the  case  where  V  =  0 

iz 

is  that  the  logarithmic  derivative  be  continuous  at  x  =*  +  d 


\Y  .  *  f  Y  /i|>  )  , 

vin  d_  1  out  out;d+ 

where  d+  =  Limfd  +  e).  This  condition  yields  the  dispersion  relation 


tan  (Q^d)  =  Qq  even  modes 


(31a) 


Q.  cotan  ( Q . d]  =  -  0  odd  modes 
i  1  i  o 


(31b) 


The  effect  of  finite  current  channel  width  are  contained  in  the  argument  of 
the  tangent  (cotangent)  function  for  the  even  modes  (odd  modes).  We 
analyze  the  dispersion  relation  for  the  even  modes  by  first  inverting  Eq. 
(31a)  to  obtain 


Qi  »  4  [ Arc  tan  f  Qq / Q ± ) ] 
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We  can  rewrice  Eq .  (32)  as 


0.  =  4  i  tan  4'  0  /0j  )  +  nnr  !  ( "*  3 ) 

•i  d  -  ° 

where  tan  [Q^/Q.J  is  the  principal  value  of  Arc  tan  iO^/O^I  and  m  is  an 
integer.  Equation  (33)  is  the  fundamental  nonlocal  dispersion  relation. 
We  note  that  when  d  =  0,  the  dispersion  relation  0Q  =  0  from  Eq.  (31a)  does 
not  give  instability  since  the  currents  driving  the  instability  vanish. 

We  can  compare  Eq.  (33)  and  Eq.  (19)  to  understand  the  local  and 
nonlocal  limits.  We  recall  that  the  local  limit  means  that  the  current 
channel  width  is  infinite  so  that  Eq.  (19)  describes  the  mode  structure 
throughout  the  space.  In  the  local  theory,  Eq.  (19)  is  also  Fourier 
analyzed  with  respect  to  x  to  obtain 


Q;  =  k 

i  x 


2  2 

and  if  we  assume  k“  »  k~  we  regain  the  earlier  results  of  Ossakow  and 

y  „  x  °  - 

Chaturvedi  (1979),  Q.-»  0,  corresponding  to  as  ai..  Comparing  Eqs.  (33) 

1  ~r 

and  (34)  we  see  that  in  the  nonlocal  analysis,  the  right  hand  side  of  Eq. 
(33)  provides  the  inverse  scale  length,  .  The  role  played  by  the 

current  channel,  then,  is  to  form  a  wavepacket  of  width  of  the  order  of 
localized  within  the  current  channel. 


In  the  limit  |0  /0,  I  »  1,  tan  *7  0  /O, )  approaches  m/2.  This 
o  i  o  i 

happens  in  the  local  limit  when  m  =  as  can  be  seen  by  dividing  Eq.  (22) 
by  Eq.  (23)  to  obtain 


tan  lfOo/Qil 


oi  f  oi  +  iv) _ 

F  0)  -  0)  )  ( OD  -  0)_  I 


As  we  approach  the  local  limit  we  have  from  Eq.  (33) 


Ch  =  [  m  +  1/2  ) 7T /d 


Figure  4 

Plot  of  growth  rate  y  versus  wavenumber  k  in  the  collisional  limit  for 
diffuse  boundary  case.  The  profiles  used  for  the  current  and  the 
density  are  given  in  Eqs.  (53)  and  (54),  respectively.  The  parameters 
used  are  the  same  as  in  Fig.  3. 


0.4 


A 

k 


Figure  3 

Plot  of  growth  rate  y  versus  k  for  d  *0.1,  1.0,  and  10.0  for  sharp 

boundary  case.  The  parameters  chosen  are  / IV^/LI  =  1, 

Q./IV./LI  =  104,  and  v  ./|V./L|  *  10~4. 
id  el  a 
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0.0030 


^in 

2i 

Figure  2 

Plot  of  growth  rate  (y/0  )  versus  the  ion-neutral  collision  frequency 

A  -4  , 

fv.  /0, )  for  k  d  *  0. 3,  V,/L  *  5.0,  and  v  /Q  =  iO  .  Curve  A  shows 
v  in  i;  y  d  ei  e 

the  exact  numerical  solutions  of  Eq.  (33).  Curves  B  and  C  represent 
Eqs.  (46)  and  (42)  which  are  solutions  in  the  inertial  and  collisional 
limits  for  long  wavelength  modes.  Curves  B'  and  C'  represent  the 
simplified  analytical  expressions,  Eqs.  (52)  and  (43),  in  the  same 
limits.  Curve  D  is  the  local  solution,  ui  ,  given  by  Eq.  (38). 
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that  both  the  long  and  short  wavelength  modes  driven  by  the  current 

convective  instability  should  have  longer  growth  times  than  previously 

estimated.  Other  data  sets  have  also  shown  very  narrow  high  latitude  F- 

region  field-aligned  currents.  From  the  results  of  this  paper,  we  expect 

the  nonlocal  nature  to  have  an  important  impact. 

Recently  Bvthrow  et  al.  (1984)  using  HILAT  satellite  data  have 

measured  very  large  downward  currents  (-  100  pA/m2),  presumably  carried  by 

cold  upward  drifting  electrons,  in  the  high  latitude  topside  F-region 

ionosphere.  The  dimensions  of  these  currents  perpendicular  to  the  magnetic 

field  is  of  order  the  transverse  scale  size  (-3  km)  of  precipitation 

induced  density  structures.  Turbulent  horizontal  ion  drifts  were  also 

measured  in  conjunction  with  the  large  currents  and  transverse  density 

gradients.  Their  data  shows  that  d  -  1.9  km,  L  =  6  km,  and  V,  *  30  km/sec. 

-4  “  -3 

In  the  case  of  ionospheric  F  region,  v  j /Q.  -  10  and  v  /ft.  -  10  .  For 

6a  o  n  i 

these  parameters  Eq .  (33)  yields  the  growth  rates  as  y  -  0.0028  0^, 

0.00260; ,  and  0.00146^.0  ,  for  X  =  1,  2,  and  10  km,  respectively.  Whereas, 

A.  A. 

for  d/L  »  1,  the  local  growth  rate  is  -  0.0031  For  narrower  current 

channel  and  large  density  gradients  or  longer  wavelength  modes  the  growth 
rate  could  be  substantially  lower. 

Finally,  several  aspects  of  the  theory  of  the  effects  of  finite 

current  channel  width  on  the  current  convective  instability  have  not  been 

discussed  here.  For  very  large  currents,  the  self-consistent  role  of 

magnetic  shear  must  be  included  [Huba  and  Ossakow,  1980].  An  approximate 

value  of  the  scale  length  for  magnetic  shear  is  L  s  cB/4irJ  with  c,  the 

s 

speed  of  light,  B  the  magnetic  field,  and  J  the  current  density. 

For  J  a  10  ^A/ra^  and  B  -  0.5  G  we  find  L  -  104  km.  For  d  <<  L  where  d  is 

s  s 

the  current  channel  width  we  do  not  expect  magnetic  shear  effects  to  play 
an  important  role.  In  addition,  we  have  not  considered  self-consistent 
velocity  shear  and  its  effects  on  the  current  convective  instability  with  a 
finite  current  channel  width,  the  finite  parallel  extent  of  the  auroral 
ionosphere,  electromagnetic  effects,  or  nonlinear  aspects.  We  are 
presently  investigating  several  of  these  aforementioned  effects. 


total  integrated  currents  in  order  to  predict  the  behaviour  of  the  olasma 
at  any  point.  Second,  because  there  is  no  plane  wave  in  the  x  direction, 
an  observer  passing  through  the  current  channel  in  the  x  direction  does  not 
see  the  usual  Doppler  shift  k  V  that  would  be  predicted  by  a  local  theorv, 

where  is  the  velocity  of  the  observer.  This  point  has  also  been  noted 

by  Ganguli  et  al .  (1984)  in  connection  with  the  current  driven  ion 

cyclotron  instability.  The  two  perpendicular  directions,  x  and  y,  are  no 
longer  equivalent. 

The  wavef  unction  as  given  by  Eas.  (21)  and  (30)  shows  that  it  is 

localized  around  the  point  x  *  0  with  a  width  Lx  *  1/0^  given  by  Eq . 
(33).  The  fastest  growing  mode  (m  =  0)  has  a  single  peak  around  x  ■  0. 
Note  that  the  width  of  the  wave  packet  is  directly  proportional  to  d,  which 
implies  that  in  the  presence  of  wider  channels  the  wavefunction  is  spread 
out,  whereas  in  the  presence  of  narrow  channels  the  wavefunction  is  sharply 
peaked . 

An  important  result  is  that  in  the  long  wavelength  limit,  the  growth 

rate  scales  as  some  power  of  the  product  of  the  magnitude  of  the  drift 

2  ct 

velocity  and  the  current  channel  width  (V-d/L")  ,  a  <  1.  This  suggests 

d 

that  the  growth  rate  would  be  the  same  whether  a  weaker  current  is 

distributed  over  a  larger  channel  or  a  stronger  current  is  concentrated  in 

a  narrow  channel.  With  regard  to  the  auroral  ionosphere,  if  the  densitv 

gradient  scale  length  is  around  20  kms,  then  we  see  from  Fig.  4,  that  the 

growth  rate  of  smaller  scale  modes,  say  \  -  1  km  corresponding  to 

k„L  '<  100,  is  reduced  from  0.3  V./L  to  0.15  V./L  by  a  current  channel  of 
7  d  d  -1-1 

width  2d  =*  40  km.  The  above  growth  rate  corresponds  to  0.75  10  s 

_  2 

if  Vj  ■  50  m/s  (corresponding  to  currents  of  order  uA/m  ).  Since 

v.  **  V./L,  these  parameters  correspond  to  the  high  latitude  F  region.  We 
m  d 

also  point  out  that  the  short  wavelength  modes  (k^L  <  1)  are  not  affected 
by  the  current  channels  of  40  kms  width. 

Vickrey  et  al.  (1980)  have  computed  the  growth  rates  from  the  current 
convective  instability  to  explain  small  scale  plasma  density  fluctuations 
near  high  latitude  F  region  large  scale  plasma  enhancements.  However, 
their  data  indicate  that  the  transverse  dimensions  of  the  currents  driving 
the  instability  are  comparable  to  the  perpendicular  (to  the  geomagnetic 
field)  scale  sizes  of  the  plasma  enhancements.  Our  results  suggest  that 
their  growth  rates  are  overestimates.  From  our  analysis  we  would  predict 


V. 


DISCUSSION  AND  SUMMARY 


We  have  presented  an  analysis  of  the  effects  of  finite  current  channel 
width  on  the  current  convective  instability.  The  current  convective 
instability  results  from  the  coupling  of  a  magnetic  field  aligned  current 
and  a  density  gradient  perpendicular  to  the  magnetic  field  in  the  presence 
of  electron-neutral  collisions.  For  a  current  with  perpendicular  scale 

size  2d  and  a  plasma  density  gradient  scale  length  L,  our  analytical 

results  indicate  a  monotonic  decrease  in  the  growth  rate  y  of  the  current 
convective  instability  for  d/L  <  1  both  in  the  collisional  fv  »  w )  and 
inertial  limits  { v .  <<  w).  For  d/L  «  1,  we  find  y  «  V,d/L^  in  the 

in  *  2  o  d 

collisional  limit  while  y  <*  fV,d/L  }  ,  a  <  1  in  the  inertial  limit.  In 

a 

general,  the  growth  rate  of  the  current  convective  instability  is  reduced 
by  the  finite  width  of  the  current  channel  from  that  of  the  local  case  in 
which  the  current  distribution  is  infinite  in  extent. 

It  is  significant  to  note  that,  in  the  long  wavelength  limit  V  «  i 
where  the  nonlocal  effects  are  most  prominent,  the  dependence  on  Vj(x) 
enters  the  nonlocal  results  in  the  sharp-boundary  case  in  the  form 
of  ( V^d) ,  the  "area"  under  the  function  V^(x)  in  the  current  channel.  It 

is  clear  that  this  quantity  is  proportional  to  the  total  current  in  the 
channel  per  unit  length  in  the  y-direction  (divided  by  the  density  of  the 
current  carriers  and  the  electric  charge).  This  point  is  to  be  contrasted 
with  the  local  dispersion  relation  (34)  in  which  the  growth  rate  depends  on 
the  current  velocity  itself,  a  local  quantity.  In  addition,  the  nonlocal 
growth  rate  depends  on  the  fractional  plasma  density  gradient  across  the 
current  channel  (i.e.  d/L).  The  growth  rate  is  reduced  because  in  a  growth 
period  (V^/L)  )  the  fluid  element  does  not  sample  the  entire  gradient 

(I/L)  but  a  smaller  fraction  of  the  density  gradient  (f/L  ,  f  »  d/L). 
Furthermore,  a  single  value  of  y  corresponds  to  the  entire  wavepacket  which 
occupies  the  current  channel.  These  nonlocal  features  are  well-known  and 
general  to  all  intrinsically  nonlocal  perturbations  such  as  the  tearing 
mode  (see,  for  example,  Chen  and  Palmadesso,  1984).  With  respect  to 
perturbations  that  have  valid  local  limits,  certain  differences  in  the 
interpretation  of  theoretical  and  observational  results  must  be  kept  in 
mind  as  one  goes  from  the  local  to  nonlocal  regime.  First,  with  finite 
current  channels,  observations  must  provide  global  quantities  such  as  the 


50".  Furthermore,  we  notice  that  the  longer  wavelength  modes  are  more 
easily  stabilized.  That  is,  a  wider  channel  (d  >  10)  significantly  reduces 
the  growth  rate  of  modes  with  wavelength  of  (k  =  0.5),  whereas  channel 
width  of  d  *  0.1  reduces  the  growth  rate  of  modes  with  wavelengths 

of  0.2nL(k  =  10)  by  the  same  amount.  From  figures  3  and  4,  we  can  conclude 
that  the  growth  rate  of  the  current  convective  instability  is  reduced  if 
the  current  channel  width  is  such  that  kd  <  1.  Similar  results  were 

"w 

noticed  in  the  context  of  collisionless  current  driven  ion-cyclotron  waves 

by  Bakshi  et  al.  (1983).  We  note  that  we  find  the  real  part  of  the 

-4 

frequency  to  be  much  smaller  than  the  growth  rate,  «-  10  <<  y. 

We  have  also  numerically  solved  equation  (11)  for  a  more  realistic 

diffuse  profile  of  density  and  current.  We  choose  a  density  profile  given 

by 


n(x) 


n 

o 


1  +  e  tanh  (x/L) 
1  -  e 


(53) 


so  that 


n*  _  j  e  sech^(x/L) 
n  L  1  +  e  tanh  (x/L) 

We  take  e  »»  0.8  so  that  n"/n  is  a  maximum  at  x/L  =  x  =»  -  0.55.  We  find 

o 

then  that  (n'/n)mo„  =  1/L.  In  addition  we  consider  a  smooth  profile  for  the 
max 

current  velocity  V,(x)z  with 

U 

V^(x)  *  exp  [-  (x-xQ)^/d^j  (54) 

We  take  the  same  parameters  as  in  the  sharp  boundary  case.  In  Fig.  5  we 
plot  y  versus  k  L  for  several  values  of  d/L.  Similar  nonlocal  behavior  is 

/  A  *i 

seen  for  the  growth  rate  y.  In  Fig.  6,  we  plot  y  versus  d/L.  The  results 
shown  in  these  figures  are  similar  to  those  obtained  by  using  a  simple 
waterbag  model  (Figs.  3  and  4)  leading  us  to  believe  that  the  simple 
waterbag  model  adequately  models  the  essential  physics,  namely  that  the 
finite  current  channel  width  has  a  stabilizing  influence  on  the  current 
convective  instability. 
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9. 


is  chac  for  d  <<  A  the  growth  rate  scales  as  [V,d/L  ,a,  where  a  <  1.  This 

a  - 

result  shows  that  in  the  long  wavelength  limit  the  growth  rate  depends  on 

the  total  current  in  the  region  where  the  wave  packet  is  localized  rather 

than  on  as  is  true  in  the  local  case. 

In  solving  equation  (33)  ,  we  consider  the  fastest  growing  mode  (m  =»  0) 

and  use  the  following  parameter  values:  v.  /  IV’./LI  *  1,  ft./|V\/Ll  =  10^ , 

^  in  d  id 

Vei^e  =  ^  and  *  ^VeiVin^e^i^  (Chaturvedi  and  Ossakow,  1981; 

Satyanarayana  and  Ossakow,  1984).  In  Fig.  3,  we  plot  the  normalized  growth 

rate  y  =  y/|V,/L|  as  a  function  of  k  =  k  L  for  various  values  of  the  width 

d  y 

of  the  current  channel  normalized  to  density  gradient  scale  length  d  =  d/L 
*  0.1,  1.0,  and  10.0.  First,  we  see  the  nonlocal  behaviour  of  the 
instability  where  the  growth  rate  is  reduced  for  k  =  k^L  <  1  and  slowly 
increases  and  saturates  to  the  local  limit  at  k  »  1  (Huba,  1984).  Second, 
the  larger  the  value  of  d,  the  smaller  is  the  k  value  for  which  the  local 

*  A  A  A  A 

limit  y  is  achieved;  for  d  =»  10,  y  approaches  yT  for  k  >  2  whereas 

A  A  A  A  L  ** 

for  d  =  0.1,  y  approaches  y  for  k  ■»  100(not  shown  in  figure).  The  local 

A  Ll 

limit,  y  ,  is  given  by  the  imaginary  part  of  Eq.  (38)  normalized  to 
Li 

I V^/L | .  Thus  we  can  conclude  that  the  current  channel  width  has  be  to 
small  (d  <  1)  to  substantially  reduce  the  growth  rate  of  the  short 

wavelength  modes  (k  >  1). 

A  A  A 

In  Figure  4  we  plot  y  as  a  function  of  d  for  k  *  0. 5,  1.0,  10.0,  and 
100.0.  As  we  discussed  earlier,  we  see  that  as  d  *  ®  we  recover  the 
conventional  current  convective  growth  rate  (y  given  in  Eq.  (38)).  The 

A  L 

second  point  is  that,  as  d  is  decreased,  the  growth  rate  decreases, 
eventually  going  to  zero  for  d  +  0.  The  quenching  of  the  instability 
should  be  expected,  because  as  d  +  0  the  currents  that  are  driving  the 
current  convective  instability  become  highly  localized  and  thus  do  not 
sample  the  entire  density  gradient.  Since  the  instability  is  driven  by  the 
density  gradient  in  the  presence  of  field  aligned  currents,  as  the  current 
channel  width  (d)  is  reduced  in  relation  to  the  density  gradient  scale 
length  (L),  the  growth  rate  of  the  instability  is  reduced.  The  instability 
is  most  sensitive  to  the  finite  width  of  the  current  channel  when  the 
current  channel  width,  d,  is  comparable  to  or  less  than  the  gradient  scale 
length,  L.  In  fact,  when  the  wavelength  is  comparable  to  the  current 
channel  width,  A  '  d  implying  kd  -  2ir ,  the  growth  rate  drops  by  at  least 


In  this  Unit  where  !Q  /Q .  |  «  1,  v.  «  ja>|  and  k“,  v.  «  |oj|,  the  growth 

'oi  in  zl  in 

rate  depends  on  the  square  root  of  the  product  of  V  and  d.  An  important 

d 

conclusion  is  that  in  both  limits  discussed  above  the  growth  rate  scales 
as  (V^d/L")01,  a  <  1.  We  note  that  the  growth  rates  given  in  Eqs.  (42)-(4A) 
can  be  shown  to  be  less  than  the  local  growth  rate,  |u)+|  ,  as  assumed 
earlier.  In  the  next  section  we  solve  numerically  the  full  nonlocal 
dispersion  relation  [Eq.  (33)]  in  order  to  study  the  growth  rate  over  a 

larger  range  of  kyL  and  d/L. 

IV.  NUMERICAL  RESULTS 

In  order  to  illustrate  the  properties  of  Eq.  (33)  and  the  scaling 

properties  of  the  growth  rate  in  the  collisional  and  inertial  domains,  we 

have  chosen,  for  convenience,  the  case  where  kzl  =  1.  This  choice 

corresponds  to  k  =  k  (v  /O  )(v.  /fi.)  «  k  .  With  this  choice,  the  local 
z  y  ei  e  in  i  y 

dispersion  relation  (equation  (38))  exhibits  a  maximum  at 
2/3  2  2  2/3 

v.  /Q.  =  (1/24)  (u)  /n,  )  .  This  is  shown  in  Fig.  2  where  we 

in  l  oi 

plot  versus  v^/fL.  As  a  numerical  example,  we  choose 

=  10  \  | V^/L|  *  5,  d/L  =  0.3  and  k^L  =  1  (Bythrow  et  al.  ,  1984). 

Curve  A  shows  the  exact  solution  of  Eq.  (33)  obtained  by  solving  it 

numerically.  Curves  B  and  C  represent  the  growth  rates  given  by  the 

analytical  expressions  in  Eqs.  (46)  and  Eq.  (42)  in  the  inertial  and 

collisional  domains,  respectively.  These  curves  show  that  there  is  good 

agreement  between  the  numerical  results  and  the  analytical  expressions 
given  in  Eqs.  (42)  and  (46).  Curves  B'  and  C'  show  the  simplified 

analytical  expressions  given  by  Eqs.  (43)  and  Eq.  (52).  Curve  D  gives  the 
local  growth  rate,  w  ,  given  by  Eq.  (38).  The  figure  shows  that  for  small 
current  channel  widths,  d  0,3  X/2tt,  where  \  is  the  wavelength  of  the 

perturbation,  the  growth  rate  is  reduced  by  an  order  of  magnitude.  For  the 

particular  choice  of  the  parallel  wavenumber  the  growth  rate  peaks  at 

-3  -3 

■j .  /n.  »  2.0  10  whereas  it  is  maximized  at  v.  /V. .  -  5.0  10  in  the 

in  1  in  i 

nonlocal  limit.  In  addition,  curves  B'  and  C'  show  that  the  Eqs.  (43) 

and  (52),  which  contain  the  scaling  properties  of  the  growth  rate,  agree 

fairly  well  with  the  exact  numerical  results 

in  and  «  0^  regimes,  respectively.  Based  on  this  agreement 

we  can  say  that  the  growth  rate  scales  as  [v.  /o,)*^  in  the  collisional 

-1/2  ln  1 

domain  and  as  (v^  /f2 ^ )  in  the  inertial  limit.  A  more  important  result 


Inertial  limit 

In  the  inertial  limit  where  vin  «  | :o t •  the  solutions  of  Eq.  (41)  are 
to ^  =  -(  iv^/2)  {l  +  lll  +  4ai“(l+<7)/v“  j^^f/fl+i^l  (46) 

_  *2  9 

where  <„  =  fl/k  d)  and  v,  =  v  +  v,  k  ,<1/2.  To  obtain  a  simple  expression 
2  v  y  ’  <  in  zl  2 

for  the  growth  rate,  we  consider  small  channel  widths,  (d  <<  L) ,  where  the 
growth  rate  is  much  smaller  than  the  local  growth  rate.  Thus  for 
m  <<  Co  Eq.  (40)  yields  in  the  inertial  limit 

a)  to  kd  =  -  to2  f  1  +  ik2,v.  /to 1 1/2  (47) 

+  -  *-  zl  in 

“2 

Equation  (47)  can  be  simplified  for  the  case  k  ^  vin  »  |to|,  i.e.  , 

fk2/k2)fn  /v  .  )  »  |  CD  |  /Q. .  The  growth  rate  can  then  be  written  as 
v  z  ’ v  e  ei '  1 

y  -  lvai/a.)l,3^V  )2/3(  <*»> 

From  Chaturvedi  and  Ossakow,  (1981)  we  have,  in  the  local  inertial  limit 
the  maximum  growth  rate 

Yt  -  !V°l/ae)‘/3  ('V/L,2/3  (4,> 

As  a  result,  we  have  for  the  nonlocal  long  wavelength  inertial  limit, 

y  *  (kyd)4’/,3rL  (50> 

Equation  (50)  shows  that  the  nonlocal  growth  rate  in  the  inertial 
limit  y  is  less  than  the  local  growth  rate  since  k^d  «  1  in  the  long 
wavelength  limit.  For  «  lwl»  i.e.,  fk2/k2)  fOe/ve^)  «  Icol/O^,  Eq. 

(47)  yields 

y  ■  (o)+to_kd)  ^2  (51) 

By  using  Eqs.  (24)  -  (27),  this  expression  can  be  rewritten  as 

T-  (k„k„/k)l/2  (Ml/2  |VL|1/2(d/L)1/2  (52) 


Collisions!  limit 


In  this  limit  where  »  |oi|  ,  Sq.(41)  has  the  following  roots: 
u^  =  -(iv>/2 }  { 1  +  [  1  +  4  u*  (  1  +  Vv^jl/?  \!\  1  +  <1)  (42) 

-o  i/2 

where  k  ■  (1  +  k“  )  /k  d  and  vv  =  v  +  v,  A  simple  expression  for 

1  zl  y  >  in  1 

the  growth  rate  can  be  obtained  by  examining  the  small  channel  width  limit, 
d  «  L,  where  the  growth  rate  is  much  less  than  the  local  growth  rate. 

Thus  for  a)  <<  o)+,  Eq.  (40)  or  (41)  yields 

*  *  *2  1/9 

m  u)  kd  *  -  icov,  f  1  +  k  ) 

+  -  in  zl 

from  which  the  growth  rate  (id  *  +  iy  ]  is  given  as 

y  =  (  kzky/k)  (1  +  kzl )_1/2  t«±/vln)  |Vd/L|  (d/L)  (43) 

It  is  worthwhile  to  compare  the  growth  rate  in  nonlocal  long  wavelength 
limit  as  given  by  Eq.  (43)  with  the  local  result  as  given  by  Chaturvedi  and 
Ossakow,  (1981) 

rL  -  (VV  <Ql/viJ  (VL|  ' 1  +  Kir'  <44) 

We  then  find  the  simple  result 

Y  =  kyd(l  +  k^)1/2YL  (45) 

Since  kyd  <  1  in  the  nonlocal  long  wavelength  limit,  the  nonlocal  growth 

rate  is  less  than  the  local  growth  rate.  Several  important  conclusions  can 

be  drawn  from  this  nonlocal  analysis.  In  this  collisional  long  wavelength 

limit,  (i)  the  growth  rate  depends  on  the  product  of  the  magnitude  of  the 

drift  velocity,  V , ,  and  the  current  channel  width,  d;  (ii)  the  growth  rate 
a 

scales  inversely  with  v^n»  the  ion-neutral  collision  frequency;  and  (iii), 
the  growth  rate  scales  inversely  with  the  gradient  scale  length,  L. 


We  set  m  =  0  and  consider  the  fastest  growing;  mode, 
definition  of  Q  from  Eq.  (23),  Eq.  (35)  becomes 


(1  -  'J)+/U)l(l  -  (D_/w) 
k  (1^ 


Then,  with  the 


(36) 


Thus,  as  we  approach  the  local  limit,  we  have  Lx  *  2d  h  and  1^  plays  the 
role  of  the  wavelength  in  the  x  direction.  For  v» ! oj I  «  1  the  new  roots 
are 


-  (iv/2)  [  1  +  (  1  +  4  J  (1  +  ir2/4k2d2)/v2  )1/2] 


ry  A  A  A  A 

(  1  +  IT  Mk/V' 


(37) 


Furthermore,  we  regain  the  local  growth  rates  for  k^d  ♦  00 ,  namely. 


r  oj+  *  (-iv/2){  1  —  |  l  +  4(a)“/v  )J  ^  1 


(38) 


In  the  limit  | Qq /Q^ I  «  1»  the  inverse  tangent  function  is 


approximately  equal  to  Q^/Q^.  Then  Eq.  (33)  becomes 


VQi 31  Qi d 


(39) 


Using  Eqs.  (22)  and  (23)  we  obtain 

*9  J  M  A  A 

(w  -  u>+)(u )  -  u_)  =  -  +  ivlnHl  +  ikz1vin/(“  +  ivin) J  /kd 

where  k2,  =  fk2/k2)fn  /v  .  )fn./v.  U  The  above  equation  is  valid  for 
„  .  zl  v  z  «  «  e  ei ' v  i  in 

f  Q  /Q  )  «  1  or  IQ.dl  «  1  which  is  essentially  the  long  wavelength  limit, 
v  o  i '  i 

kyd  «  1.  Equation  (40)  can  then  be  solved  both  in  the  highly  collisional 
and  weakly  collisional  domains.  Using  Eqs.  (24)  and  (25),  Eq.  (40)  can  be 
written  as 

U)2  +  iv  m  +  u)2  *  +  o)(u>  +  iv  ^ n H  1  +  *  kzlvin^  w  +  ^vin^  /kd  (^D 


The  negative  (positive)  sign  on  the  right  hand  side  yields  a  growing 
(damped)  mode.  In  the  following  anlaysis  we  consider  the  growing  modes. 
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LOS  ALAMOS,  NM  85544 

OICY  ATTN  DOC  CON  FOR  J.  RREEDLO'YE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  TECH  INFO 
DEPT 

OICY  ATTN  DOC  CON  FOR  L-389  R.  OTT 
OICY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 


OICY 

ATTN 

DOC 

CON 

FOR 

J  . 

WOLCOTT 

OICY 

ATTN 

DOC 

CON 

FOR 

R. 

F.  TASCHEK 

OICY 

ATTN 

DOC 

CON 

FOR 

E. 

JONES 

OICY 

ATTN 

DOC 

CON 

FOR 

J  . 

MALIK 

OICY 

ATTN 

DOC 

CON 

FOR 

R. 

JEFFRIES 

OICY 

ATTN 

DOC 

CON 

FOR 

J  . 

ZINN 

OICY 

ATTN 

DOC 

CON 

FOR 

P. 

KEATON 

OICY 

ATTN 

DOC 

CON 

FOR 

D. 

WESTERVELT 

OICY 

ATTN 

D. 

SAPPFNFIELD 

S  AN  0 1  A  LABORATORIES 
P.O.  BOX  5800 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  W.  RROWN 
OICY  ATTN  DOC  CON  FOR  A. 

THORNBROUGH 

01C Y  ATTN  DOC  CON  FOR  T.  WRIGHT 

0 1 C  Y  ATTN  DOC  CON  FOR  D.  DAHLGREN 

OICY  ATTN  DOC  CON  FOR  3141 

OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT 
DIV 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  B.  MURPHEY 

OICY  ATTN  DOC  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  DR.  YO  SONG 


OTHE  R  GOVERNMENT 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMNI CATIONS  &  INFO 
ADMIN 

30ULDER,  CO  30303 


01  CY 

ATTN 

A. 

JEAN  (UNCLASS  ONLY) 

01  CY 

ATTN 

W. 

UTLAUT 

01  CY 

ATTN 

D. 

CROMBIE 

01  CY 

ATTN 

L. 

BERRY 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
535  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
01CY  ATTN  D.B.  COX 
01  C  Y  ATTM  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
01 C  Y  ATTN  G.  HYDE 


NATIONAL  OCEANIC  S  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOUuDER,  CO  80302 
01CY  ATTN  R.  GRU88 
01CY  ATTN  AERONOMY  LAS  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 


AEROSPACE  CORPORATION 
P.O.  BOX  92957 


LOS  ANGELES, 

CA 

90009 

01  CY 

ATTN 

I. 

GARFUNKEL 

01  CY 

ATTN 

T. 

SALMI 

01  CY 

ATTN 

V. 

JOSEPHSON 

01  CY 

ATTN 

S. 

BOUER 

01  CY 

ATTN 

D. 

OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
8URLINGT0N,  MA  01303 

01CY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 


CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 

ITHACA,  NY  14850 

01CY  ATTN  D.T.  FARLEY,  JR. 

EL£CTROS°ACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
01CY  ATTN  H.  LOGSTON 
01CY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  W i l shire  Blvd. 

Santa  Monica,  Calif  90401 
01 CY  ATTN  C.8.  GABBARD 
01 CY  ATTN  R.  LELEVIER 

ESL,  INC. 

495  JAVA  DRIVE 

SUNNYVALE,  CA  94086 
0 1 C Y  ATTN  J.  ROBERTS 
01 CY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 

SPACE  DIVISION 

VALLEY  FORGE  SPACE  CENTER 


AUSTIN,  TX  73758 


01  CY 

ATTN 

L. 

SLOAN 

01  CY 

ATTN 

R. 

THOMPSON 

8ERKELEY 

RESEARCH  ASSOCIATES 

P.O.  BOX 

983 

BERKELEY 

,  CA 

94701 

01  CY 

ATTN 

J  . 

WORKMAN 

01  CY 

ATTN 

C. 

PRETTIE 

01  CY 

ATTN 

s. 

8RECHT 

BOEING  COMPANY, 

THE 

P.O.  BOX 

3707 

SEATTLE, 

WA 

931  24 

01  CY 

ATTN 

G. 

KEISTER 

01  CY 

ATTN 

D. 

MURRAY 

01  CY 

ATTN 

G  • 

HALL 

01  CY 

ATTN 

J  . 

KENNEY 

GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 
0 1 C Y  ATTN  M.H.  BORTNER 
SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 

01CY  ATTN  F.  RE  I  BE  RT 

GENERAL  ELECTRIC  TECH  SERVICES 
CO.,  INC. 

HMES 

COURT  STREET 
SYRACUSE,  NY  13201 

0 1 C  Y  ATTN  G.  MILLMAN 


GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 

01 CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
0 1 C  Y  ATTN  TECHNICAL  LIBRARY 
0 1 C Y  ATTN  NEAL  BROUN  (UNCLASS  ONLY) 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

0 1 C Y  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01 CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

0 1 C Y  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
0  1  C Y  ATTN  J.M.  AEIN 
0 1 C  Y  ATTN  ERNEST  BAUER 
01 CY  ATTN  HANS  UOLFARD 
0 1  C Y  ATTN  JOEL  BENGSTON 

INTL  TEL  S  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

01 CY  *ATTN  TECHNICAL  LIBRARY 

J  AYCOR 

11011  TORREYANA  ROAD 

P.O.  BOX  85154 

SAN  DIEGO,  CA  92138 

0 1 CY  ATTN  J.L.  SPERLING 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

0 1  C Y  ATTN  DOCUMENT  LIBRARIAN 
0 1 CY  ATTN  THOMAS  POTEMRA 
0 1 C  Y  ATTN  JOHN  DASSOULAS 


KAMAN  SCIENCES  CORP 
P.O.  POX  7463 

COLORADO  SPRINGS,  CO  30933 
0 1 C  Y  ATTN  T.  MEAGHER 

KAMAN  T£M°0-CENTfR  FOR  ADVANCED 
STUDIES 

816  STATE  STREET  (P.O  DRAWER  00) 

SANTA  BARBARA,  CA  9310? 

01  CY  ATTN  DAS  I  AC 
01CY  ATTN  UARREN  S.  KNAPP 
01CY  ATTN  WILLIAM  MCNAMARA 
01CY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 

0 1 C  Y  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  8  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
01CY  ATTN  DEPT  60-12 

0 1 C  Y  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  %  SPACE  CO.,  INC. 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

01 CY  ATTN  MARTIN  WALT  DEPT  52-12 
01CY  ATTN  W.L.  IMHOF  DEPT  52-12 

01CY  ATTN  RICHARD  G.  JOHNSON 

DEPT  52-12 

01CY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 

01CY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

01 CY  ATTN  DAVID  M.  TOULE 
01 C  Y  ATTN  L.  LOUGHLIN 
01CY  ATTN  D.  CLARK 


T 


r 


MCDCNNEL  DOUGLAS  CORPORATION 
5301  80LSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CT  ATTN  N.  HARRIS 
01CY  ATTN  J.  MOULE 
01CY  ATTN  GEORGE  MROZ 
01CY  ATTN  W.  OLSON 
0 1 C Y  ATTN  R.U.  HALPRIN 
01 CY  ATTN  TECHNICAL 

LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 

735  STATE  STREET 

SANTA  BARBARA,  CA  93101 


Cl  CY 

ATTN 

P.  FISCHER 

01  CY 

ATTN 

W.F.  CREVIER 

01  CY 

ATTN 

STEVEN  L.  GUTSCHE 

01  CY 

ATTN 

R.  BOGUSCH 

01  CY 

ATTN 

R.  HENDRICK 

01  CY 

ATTN 

RALPH  KILB 

01  CY 

ATTN 

DAVE  SOWLE 

01  CY 

ATTN 

F.  FAJEN 

01  CY 

ATTN 

M.  SCHEIBE 

01  CY 

ATTN 

CONRAD  L.  LONGMIRE 

01  CY 

ATTN 

B.  WHITE 

01  CY 

ATTN 

R.  STAGAT 

MISSION 

RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 
ALBUQUERQUE,  NEW  MEXICO  87106 
01CY  R.  STELLINGWERF 
01CY  M.  ALME 
01 CY  L.  WRIGHT 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

0 1 C  Y  ATTN  JOHN  MORGANSTERN 
01 CY  ATTN  G.  HARDING 
0 1 C  Y  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1320  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
0 1 C  Y  ATTN  W.  HALL 
0 1 C  Y  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

0 1  C  Y  ATTN  E.C.  FIELD,  JR. 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
( NO  CLASS  TO  THIS  ADDRESS) 

01 CY  ATTN  IONOSPHERIC  RESEARCH  LAP 

PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

0 1 C Y  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

0 1 CY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  »  D  ASSOCIATES 

P.O.  SOX  9695 

MARINA  DEL  REY,  CA  90291 


01  CY 

ATTN 

FORREST  GILMORE 

0 1  C  Y 

ATTN 

WILLIAM  9.  WRIGHT 

01  CY 

ATTN 

WILLIAM  J.  KARZAS 

01  CY 

ATTN 

H.  ORY 

01  CY 

ATTN 

C.  MACDONALD 

01  CY 

ATTN 

R.  TURCO 

01CY 

ATTN 

L.  DeRAND 

01  C  Y 

ATTN 

W.  TSAI 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01 CY  ATTN  CULLEN  CRAIN 
0 1 C Y  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

0 1 C  Y  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
330  WEST  42nd  STREET 
NEW  YORK,  NY  10036 

01 CY  ATTN  VINCE  TRAPANI 
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SCIENCE  APPLICATIONS,  INC. 

1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

OUY  ATTN  LEWIS  M.  LINSGN 
01CT  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  E.  FRIEMAN 
01CY  ATTN  E . A .  STRAKER 
01CY  ATTN  CURTIS  A.  SMITH 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
0  1  C  Y  J.  COCKAYNE 
0 1 C  Y  E.  HYMAN 

SRI  INTERNATIONAL 
333  RAVENS WOOD  AVENUE 
MENLO  PARK,  CA  94025 
0 1 C  Y  ATTN  J.  CASPER 
01CY  ATTN  DONALD  NEILSON 
01CY  ATTN  ALAN  BURNS 
01 CY  ATTN  G.  SMITH 
01CY  ATTN  R.  TSUNODA 
0 1 CY  ATTN  DAVID  A.  JOHNSON 
01CY  ATTN  WALTER  G.  CHESNUT 
OUY  ATTN  CHARLES  L.  RINO 
OUY  ATTN  WALTER  JAYE 
0 1 CY  ATTN  J.  VICKREY 
OUY  ATTN  RAY  L.  LEAD  ABR  AND 
0 1 CY  ATTN  G.  CARPENTER 
01CY  ATTN  G.  PRICE 
0 1 CY  ATTN  R.  LIVINGSTON 
01CY  ATTN  V.  GONZALES 
01CY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

01CY  ATTN  W.P.  80QUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
01CY  ATTN  JOHN  ISE,  JR. 

0 1 C  Y  ATTN  JOEL  GAR8ARIN0 

TRW  DEFENSE  3  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  8EACH,  CA  90278 
Q1CY  ATTN  R.  K.  PLEBUCH 
01CY  ATTN  S.  ALTSCHULER 
OUY  ATTN  D.  DEE 
OUY  ATTN  D/  STOCKWELL 
SNTF/1575 


VISIDYNE 

SOUTH  3EDF0RD  STREET 
BURLINGTON,  MASS  01803 
01CY  ATTN  W.  REIDY 
OUY  ATTN  J.  CARPENTER 
01CY  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

01CY  ATTN:  N.  ZABUSKY 


IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 


PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF 
NOTED) 

NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 
Code  4100 

Dr.  H.  GURSKY  -  CODE  4100 
Dr.  P.  GOODMAN  -  CODE  4180 

A.F.  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  MA  01731 
DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J  .M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  U.  BURKE 
DR.  H.  CARLSON 
DR.  J.  JASPERS 
Dr.  F.J.  RICH 
DR.  N.  MAYNARD 

BOSTON  UNIVERSITY 
DEPARTMENT  OF  ASTRONOMY 
BOSTON,  MA  02215 
DR.  J.  AARONS 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  W.E.  SWARTZ 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
OR.  M.8.  McELROY 
DR.  R.  LINDZEN 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY /NAVY  DRIVE 
ARLINGTON,  VA  22202 
OR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF 
TECHNOLOGY 

PLASMA  FUSION  CENTER 
LIBRARY,  NW16-262 
CAMBRIDGE,  MA  02139 


THE  FOLLOWING  PEOPLE  (UNLESS  OTHERWISE 


NASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREEN8ELT,  MD  20771 
DR.  K.  MAEDA 
DR.  S.  CURTIS 
DR.  M.  DU8IN 

COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUSA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO,  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND 

ENVIRONMENTAL  LABORATORY 
BOULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  G.  JOINER 

PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  °ARK,  PA  16802 
DR.  J.S.  NIS8ET 
DR.  P.R.  ROHRBAUGH 
DR.  L.A .  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  F.  KLEVANS 

SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
DR.  D.A.  HAMLIN 
DR.  E.  FRIEMAN 


STANFORD  UNIVERSITY 
STANFORD,  CA  34305 
DR.  P.M.  BANKS 

U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABEROEEN,  MD 

DR.  J.  HEIMERL 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 
DR.  L.E.  LEE 

UNIVERSITY  OF  CALIFORNIA, 
BERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 

UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-10,  MS-664 
LOS  ALAMOS,  NM  87545 
DR.  M.  PONGRATZ 
DR.  D.  SIMONS 
DR.  G.  3ARASCH 
DR.  L.  DUNCAN 
DR.  P.  8ERNHARDT 
DR.  S.P.  GARY 

UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 

JOHNS  HOPKINS  UNIVERSITY 
APPLIEO  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 
DR.  R.  GREENUALD 
DR.  C.  MENG 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 
DR.  N.  ZA8USKY 
DR.  M.  BIONDI 
DR.  E.  OVERMAN 


UNIVERSITY  OF  TEXAS 
AT  DALLAS 

CENTER  FOR  RESEARCH  SCIENCES 
P.O.  BOX  688 
RICHARDSON,  TX  75080 
DR.  R.  HEELIS 
DR.  W.  HANSON 
dr.  j.p.  McClure 

UTAH  STATE  UNIVERSITY 
4TH  AND  8TH  STREETS 
LOGAN,  UTAH  84322 
DR.  R.  HARRIS 
DR.  K.  BAKER 
DR.  R.  SCHUNK 
DR.  J.  ST. -MAURICE 

PHYSICAL  RESEARCH  LABORATORY 
PLASMA  PHYSICS  PROGRAMME 
AHMEDABAD  380  009 
INDIA 

P.J.  PATHAK,  LIBRARIAN 

LABORATORY  FOR  PLASMA  AND 
FUSION  ENERGY  STUDIES 
UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20742 
JHAN  VARYAN  HELLMAN, 
REFERENCE  LIBRARIAN 


END 


FILMED 

4-85 
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